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Perturbative QCD analysis of the nucleon's Pauli form factor F 2 (Q 2 ) 

A.V. Belitsky, Xiangdong Ji, and Feng Yuan 
Department of Physics, University of Maryland, MD 201^2-^111, College Park, USA 

We perform a perturbative QCD analysis of the nucleon's Pauli form factor F2(Q 2 ) in the asymp- 
totically large Q 2 limit. We find that the leading contribution to F2(Q 2 ) has a l/Q 6 power behavior, 
consistent with the well-known result in the literature. Its coefficient depends on the leading- and 
subleading- twist light-cone wave functions of the nucleon, the latter describing the quarks with 
one unit of orbital angular momentum. We also derive at the logarithmic accurary the asymptotic 
scaling F2(Q 2 )/ Fi(Q 2 ) ~ (log 2 Q 2 / A 2 )/Q 2 which describes recent Jefferson Lab data well. 

PACS numbers: 13.40.Gp, 12.38.Bx 
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The electromagnetic form factors are fundamental ob- 
servables of the nucleon containing important informa- 
tion about its internal nonperturbative structure. Since 
the first measurement in the mid 1950s, experimental 
studies of these observables have become an active fron- 
tier in nuclear and particle physics. More recently, with 
the development of novel experimental techniques, the 
magnetic form factor of the proton, and electric and mag- 
netic form factors of the neutron have been measured 
with unprecedented precision at Jefferson Lab and other 
facilities around the world |]|. In particular, the Pauli 
form factor of the proton F 2 extracted from the recoil 
polarizations at Jefferson Lab reveals a significant dif- 
ference from previous data and theoretical expectations 
0, 0] • The data shows that the ratio of Dirac to Pauli 
form factor Q 2 F 2 (Q 2 ) / Fi(Q 2 ) continues climbing at the 

12 !, 



largest Q s measured and seems to scale as \JQ 2 - The 
result has sparked myriad speculations on its implication 
about the underlying microscopic structure of the proton 
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FIG. 1: A leading QCD diagram contributing to the nucleon's 
form factors. 

A sim ple dimensional counting rule was devised in 
Refs. |jl Il0j| to determine the dominant power contri- 
bution to hadronic form factors at large momentum 
transfer. For the hadron helicity-conserving form factor 
Fi{Q 2 ), it predicts a dominant scaling behavior l/Q 4 . 
The power counting can be justified by QCD factor- 
ization theorems which separate short-distance quark- 
gluon interactions from soft hadron wave functions 0, 
El [H H G3 H [H E3, see Fig. □ Recently, a 
more sophisticated formalism has been developed to in- 
clude higher-order perturbative QCD (PQCD) resum- 
mation, or Sudakov form factor, which treats contri- 



butions from the small- a; partons more appropriately 
13 El El El E3 The procedure extends the appli- 
cability of leading-order PQCD predictions to moderate 
Q , but does not ch ang e the dominant power-law behav- 
ior significantly |17L I22L l23| . 

The literature on PQCD studies of F 2 (Q 2 ), however, 
is much meager. Since F 2 (Q 2 ) is related to the hadron 
helicity-flip amplitude, its power behavior is suppressed 
compared to Fi(Q 2 ) ~ 1/Q 4 , and a natural expectation 
is then 1/Q 6 . This is confirmed by calculations where 
the quark masses are the mechanism for the helicity flip 
0, Ej|- Since the up and down current quark masses 
are negligible, the dominant mechanism for the spin flip 
in QCD come from the quark orbital angular momen- 
tum and the polarization of an extra gluon which must 
be included to maintain gauge invariance j^. A gen- 
eralized power counting including the parton orbital an- 
gular momentum validates the expected scaling law for 
F 2 {Q 2 ) HI. An actual PQCD calculation for F 2 {Q 2 ) 
requires the extension of the usual formalism to include 
the quark orbital angular momentum, a technology has 
not yet been systematically developed in the literature. 
The recent development in classification of the light-cone 
wave functions [271 i=M provides the necessary ingredient 
to perform these calculations. 

In this paper we report on a first PQCD calculation 
of F 2 (Q 2 ) at asymptotically large momentum transfer 
Q 2 . Modulo logarithms, we confirm that the leading con- 
tribution to F 2 (Q 2 ) goes like 1/Q 6 . The coefficient of 
the leading-power term depends on leading order (twist- 
three) and next-to-leading order (twist-four) light-cone 
wave functions. The latter is the probability amplitude 
for one of the quarks to carry one unit of orbital angular 
momentum. We compute the leading-order perturbative 
kernel, and estimate the coefficient using models for the 
light-cone wave functions. We comment on the role of 
the Sudakov form factor in regulating possible end-point 
singularities in the phase space integrals. We also de- 
rive with the logarithmic accuracy the asymptotic scal- 
ing F 2 (Q 2 ) / Fi(Q 2 ) ~ (log 2 Q 2 /A 2 )/Q 2 which we find to 
describe the data unexpectedly well. 

To start with, we choose a coordinate system for the 
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process P + 7* — > P' such that P and P' form conju- 
gate light-cone vectors. We let the light-cone compo- 
nents P + and P'~ to be large in the asymptotic limit, 
Q 2 = 2P+P'~. The F 2 (Q 2 ) form factor can be extracted 
from the helicity-flip matrix element 
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where is the quark electromagnetic current and u(P) 
is an on-shell spinor of the nuclcon with momentum P. 
In the above coordinate system, the spin of the nucleon 
actually remains the same in the initial and final states. 

A perturbative analysis of the above matrix element 
is done by computing leading-order Feynman diagrams 
with two gluon exchanges, for a typical graph see Fig. ^ 
in which the initial (final) nucleon contains three quarks 
with longitudinal momenta xiP (y%P') and transverse 
momenta ki (fe£) of order Aqcd- Two hard-gluon ex- 
changes ensure that the three quarks in the final state 
propagate collinearly after the injection of a large mo- 
mentum transfer q^. Since up and down quarks are 
light, quark helicity is approximately conserved during 
the hard scattering. Therefore, to produce a nucleon-spin 
flip, the quark orbital angular momentum in the initial 
and final states must differ by one unit. The leading 
contribution to F2(Q 2 ) comes from the configurations in 
which the quarks in the initial or the final state carry zero 
unit, and those in the other state carry one unit, of or- 
bital angular momentum. For definiteness, let us assume 
that the final state quarks are symmetric in azimuth. 

To isolate the leading contribution, we expand the hard 
part of the diagram in k 2 /Q 2 in the limit of large Q 2 and 
Xi 7^ (a procedure dubbed the collinear expansion). 
We will comment on the x% — ► case later. Since the 
final state nucleon has no orbital angular momentum, we 
throw away all the subleading terms in k' t . On the other 
hand, the quarks in the initial state nucleon have one unit 
of orbital angular momentum, and hence the hard part 
must have linear terms in the quark transverse momenta. 
Therefore, the leading hard part, we are interested in, has 
a structure kiT(xi,x 2 ,x 3 ,yi,y 2 ,y3;Q 2 )- 

Let us consider the simplification of the nucleon wave 
functions after the collinear expansion. For definiteness, 
we consider the proton form factor. The light-cone wave 
function for the final state is 
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[dx][d 2 k] 



V>l(Kl, «2, K3) 



(2) 



Vl / 2 12 J ^xTx^Ei 
xe ahc u^{ Kl ) (4 i ( K2 )d t cT ( K3 ) - 4^2)^3)} |0) , 

where the argument K,^ IS Si shorthand notation for 
(xi, ki), and the integration measures for the quark mo- 
menta are 

[dx] = dxi dx2 dx 3 5 (xi + X2 + x 3 — 1) , 
[d 2 k] = d 2 ki d 2 k 2 d 2 k 3 6 (2 *> (fei +k 2 + fe 3 ) • 



If the hard part has no dependence on fe£, we can ig- 
nore the transverse momentum dependence in the quark 
creation operators and define the twist-three amplitude 

^3(^1,^2,^3) = 2 J[d 2 k] ipi(ni, K 2 , K3) • 

Then the final nucleon state can be simplified to, 

[dx] 



i^i/ 2 - 24 y 



:@ 3 (xi,X 2 ,X 3 ) 



y/XiX2X 3 

abc ul^ Xl ) {ul L (x2)d^(x 3 ) - 4^2)^3)} |0) 



(3) 



In the first quantization formalism, the second line can 
be replaced by the standard SU (6) wave function. 

Since the quarks in the initial nucleon state must have 
a total helicity 1/2, only the following wave function com- 
ponent is relevant [^Eif 
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[dx][d 2 k] 



{&1J_V>3 + k 2 ±1p4,} «2, K3) 



xe 



^/XlX 2 X 3 

abc ^ T ( K2 ){^ i ( K1 ) d t T(K3) _ d t i(Kl)u t T(K3) || 0)i(4) 



where we use the notation k± = k x — ik y . 

The spin-isospin structure of this wave function is ex- 
actly the same as that in Eq. J2J. After the collinear 
expansion, the above wave function will be convoluted 
with transverse momenta of quarks. We define the twist- 
four amplitudes via 

f [d 2 k] 

$4,{X2,Xi,X 3 ) = 2 I — — k 3 •{fcj.V's + ^2^4} «2, K3), 
J Mx 3 



W i (x u x 2 ,x 3 ) = 2 



[d 2 k] 
Mx 2 



k 2 -{fcl^3 + k 2 lf>i} K2, K 3 ). 



Apart from a gluon-potential dependent term, the above 
light-cone amplitudes are the same as those defined in 
Ref. jl^l . In the following, we include the former so that 
the result is gauge invariant. To the order that we are 
working, the gluon potential terms arise from perturba- 
tive diagrams with explicit gluons attached to the nucleon 
blobs. These contributions are generally considered as 
dynamically suppressed [30| . 

We introduce effective nucleon wave functions for the 
initial state after integrating over the transverse mo- 
menta weighted with a momentum factor from the hard 
part, 
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[dx 



:X1^ , 4{X2,X1 1 X 3 ) 



^JX\X2X 3 

xe abc u^{x x ) {4^x2)4^x3) - 4 ; (z 2 )4 T (:r3)} |0) 



(5) 



ln>i /2 [fca] = g 



[dx] 



:X 3 <P4(X1,X 2 ,X 3 ) 



y/XiX 2 X 3 

Xe abc ul r ( Xl ) {4 l (x 2 )d^(x 3 ) - dl l (x 2 )4 r (x 3 )} |0) 



(6) 



where the argument on the left-hand side indicates the 
momentum being averaged. Averaging over fc 2 does not 
yield independent information. 

Using the twist-3 and -4 amplitudes, we can write the 
F2{Q 2 ) form factor in the following factorized form, 



F- 



(Q 2 ) = l[dx][dy]{x 3 ^ i (x 1 ,x 2 ,x 3 )T^({x},{y}) (7) 

+ Xl&i(x2,Xl,X 3 )T$r({x}, {y})\@ 3 (yi, V2, 2/3) , 

where {x} = (x±, x%, x 3 ). Let us see how to extract the 
hard part, T, in Fig.^ 

We use the nomenclature and strategy of Ref. ^lj by 
calling the top quark line 1, the middle 2, etc., without 
committing them to a specific flavor. Given the spin- 
isospin structure of the initial and final state nucleon 
wave functions, we assume that the first and third quarks 
have spin up and the second one spin down. Call the 
amplitudes Tj when the electromagnetic current acts on 
particle i. Clearly because of symmetry, T3 can be ob- 
tained from T\ by a suitable exchange of variables (see 
below). Using the SU(6) wave function and the quark 
charge weighting, we find the flavor structure of the hard 
part for the proton, 

TP({*}, M) = + ^±^(T 2 + T 3 ) 



(8) 



where we have omitted the argument of Ti on the right- 
hand side, and T- has y\ and y 3 interchanged. For the 
neutron, u <-> d. 

According to the above convention, we find the contri- 
bution to T\ from Fig. ^ 



Tft{{x},{y}) = -^(4™ s ) 2 
T?i-\{x},{y}) = -^(4™ s ) 2 



f 



x\x 3 yiy 3 



(9) 



where here and everywhere x = 1 — x. Cb = 2/3 
is a color factor and the contribution from the twist-4 
wave functions of the final-state nucleon has also been in- 
cluded. The corresponding contribution to T 3 can be ob- 
tained by interchanging the labels, T3J' 1 = Tfj ^(l 3), 
T 3 $ l — T^' X {1 «-* 3). This in fact is a general feature. 

To find a complete expression for the hard part, we 
must calculate perturbative diagrams displayed in Fig. 
El A straightforward evaluation yields, 



M 2 C 2 

T w ({x},{y}) = — ^{Ana s f x^ 



Q 6 



G12), (10) 



T 19 ({x},{y}) 



M 2 r 2 

-^(47ra s ) 2 [- Xl G n - XiGn] 



m 2 c, 2 

T 2 *({x},{y}) = -^(^s 



Q 6 

x [x 3 (G 22 - G 2 i - G22) - x 3 G 2 i\ 



FIG. 2: Perturbative diagrams contributing to the hard part 
of F2. Mirror symmetric graphs have to added. 



and T 2 <g = ?2^(1 <-> 3), T 3 & — Ti#(l <-* 3), and T 3 $ 
Ti^(l <-> 3). The functions Gy are defined as 
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x x x2x\y2y\y 3 



G 



x 2 x\y 3 yl x 2 x 3 x 2 y 2 y 2 
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x\x 3 x 3 yly 3 yi 



Xxx\x2y\y2 x\x\y\y 3 y\ 



(11) 



and G21 = G?ai(l «-» 3), G 22 = G 22 (l «-» 3). 

To determine the normalization for F 2 (Q 2 ), we need to 
know the light-cone distribution amplitudes <P 3 , #4 and 
S?4, which can only be obtained by solving QCD nonper- 
turbatively. However, the scale evolution of these ampli- 
tudes selects at the asymptotically large Q 2 the leading 
component with a fixed small-Xi behavior. For example, 
the asymptotic form of <P 3 is x \X2X 3 whereas that of 
<p4 is x\X2- In Ref. [2^| a set of phenomenological ampli- 
tudes satisfying these asymptotic constraints have been 
proposed on a basis of conformal expansion. 

With the above wave functions, the integrals over mo- 
mentum fractions X4 and yt have logarithmic singulari- 
ties, indicating that the factorization breaks down when 
one of the quarks in the wave function becomes soft 
[l9L |2C| . It has been suggested that the higher-order 
PQCD resummation, or the Sudakov form factor, sup- 
presses the contribution at small- a; and provides an ef- 
fective cut-off for the integrals at x ~ A 2 /Q 2 , where 
A is a soft scale related to the size of the nucleon 
[Hill H3. The outcome is that the x 4 integra- 
tions contribute an extra Q 2 -dependent factor In 2 Q 2 /A 2 , 
compared to Fi(Q 2 ). Physically, the end-point diver- 
gencies indicate that quarks with different rapidity con- 
tribute equally to the hard scattering. Since the con- 
tribution from quarks with very large rapidity (small- 
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x) is suppressed by the Sudakov form factor, this Q 2 - 
dependence reflects simply the kinematic broadening of 
the quark (and gluon ) rapidity range with increasing nu- 
cleon momentum. 

For an estimate, we use asymptotic wave functions |29j . 

<£ 3 = UOx^xs/n, $4 = 12 XiX 2 (f N + Xi) , 
$ i = 12xiXz{f N -\ 1 ), (12) 

with f N = 5.3 ■ 1CT 3 GeV 2 and Ai = -2.7 • 10~ 2 GeV 2 . 
With a choice of A = (0.3 GeV) 2 , Q & F%(Q 2 ) is roughly 
0.6 GeV 6 for Q 2 = (5 - 20) GeV 2 , about 1/3 of the 
Jefferson Lab data at Q 2 = 5 GeV 2 01 ■ Of course, to 
get a more realistic PQCD prediction in this regime, one 
must have the quark distribution amplitudes appropriate 
at this scale. However, from the comparison between the 
data and PQCD predictions for F\(Q 2 ) ||, we believe 
that asymptotic PQCD is unlikely to be the dominant 
contribution to ^(Q 2 ) at Q 2 = 3 — 5 GeV 2 : one must 
take into account higher-order corrections and higher- 
twist effects. 
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FIG. 3: JLab data plotted in terms of the leading PQCD 
scaling. The low, middle, and upper data points correspond 
to A = 200, 300, 400, respectively. 

Coming back to the scaling behavior of the ratio 
-FMQ 2 ) / 'F\(Q 2 ) for which the Jefferson Lab data has 
stimulated much discussion in the literature. PQCD 
predicts the power-law scaling 1/Q 2 . With the new 
result for ^(Q 2 ), we can determine its scaling up to 
logarithmic accuracy. The strong coupling constant in 
the ratio simply cancels. The wave function evolution 
yields a factor of af /m {Q 2 ) for F\(Q 2 ) and aJ m (Q 2 ) 
for F 2 (Q 2 ) from the leading non- vanishing contribution, 
where (3 = 11 - 2ra//3. Thus PQCD predicts that 
(Q 2 /\n 2+8/m Q 2 /A 2 )(F 2 (Q 2 )/F 1 (Q 2 )) scales as a con- 
stant at large Q 2 , 8/(9/3) <C 1. Surprisingly, the Jefferson 
Lab data plotted this way (ignoring the small 8/(9(3)) 
exhibits little Q 2 variation for a range of choices of A 
as shown in Fig. [3J Since we do not expect the asymp- 
totic predictions for Fi^iQ 2 ) to work at these Q 2 , the 



observed consistency might be a sign of precocious scal- 
ing as a consequence of delicate cancellations in the ratio. 
A more detailed discussion on this issue along with more 
thorough phenomenological analyses will be given in a 
separate publication. 
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